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Techniques are developed which permit general purpose structural 
analysis computer programs to be used to generate the equations of motion 
necessary for limiting performance studies. The limiting performance 
characteristics of a system are useful In analyrlng the optimal behavior 
of a dynamic system. In particular, the limiting performance charac- 
teristics are the essential Ingredients In an efficient optimal design 
method for Isolation systems. 
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INTRODUCTION 


The use of conventional methods of optimal design has ao far been 
limited to simple systems with few design parameters. This is primarily 
due to the overwhelming computational burdens of the methods which 
Invariably require solving the system dynamics repeatedly. Recently a 
new approach that differs from the conventional methods In design 
methodology has beer proposed and extensively explored (Reference 1, 2). 

The new approach, called the indirect synthesis method, selects the design 
parameters on the basis of a limiting performance study of the dynamic 
system being designed. In the process the system dynamics need be solved 
only once and thus the computational effort is greatly reduced. The 
method has been successfully applied to problems ranging from an infinite 
degree of freedom system with two design parameters (a beam) to a five 
degree of freedom system with six design parameters (an automobile model) 
(Reference 2). Frequently, only the limiting performance calculations are 
made for an isolation system as this provides useful information concern- 
ing the optimal system response. In particular, the limiting performance 
characteristics can be used as a yardstick for gauging the desirability 
of candidate system designs. 

Several computer programs are now available for limiting performance 
studies. These are reviewed in Reference 3. Some of these programs, e.g., 
PERFORM (Reference 4) for transient systems and SYSLIPEC (Reference o) for 
steady state systems, are virtually unlimited in terms of size and complexity 
of acceptable dynamic systems. PERFORM and SYSLIPEC are available through 
COSMIC as programs LAR-11930 and LAR-11931, respectively. The limiting 
performance computer programs usually require mass, stiffness, damping, and 
other characteristic matrices as input. The program then computes the 
limiting performance characteristics or. the basis of this system information. 



In view of the success as well as the extensive use of the finite 
element (IT) method in analyzing complex structures, it is 
important to develop the methodology necessary to couple available 
general purpose finite element structural programs to a limiting per- 
formance capability. This is the purpose of the work reported here. 

Once these general purpose programs can be used to develop input for a 
limiting performance computer program, then the limiting performance 
and indirect synthesis can be applied to any system for which the general 
purpose program FE is appropriate. 

Approaches for coupling general purpose programs to limiting per- 
formance capabilities will be considered. Primary emphasis is given to 
the use of the general purpose program to develop equations of motion in 
a form that can be used by the limiting performance program. 

BACKGROUND 

Consider a linear dynamic system whose equations of motion can be 
written in matrix form as 

[M] {x} + [C] {x} + (KHx) - [F] { f } (1) 

where 

{x} ■ displacement vector 

[M] ■ NxN Mass matrix, N being the number of degrees of freedom (DOF) 
[C] “NxN damping matrix 

[K] “NxN stiffness matrix 

[F] - NxL coefficient matrix associated with the forcing function 
vector { f } , L being the number of forcing functions. 

The optimal design problem is to choose portions of the system so that 

some index of performance is minimzed (or maximized) and certain constraints 
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are satisfied. It is assumed that both the performance index (tj/) and the 


constraint functions C^(x,t) on the response variables Involve 
peak response variables, e.g., maximum stresses, displacements , ac- 
celerations. The performance index and the constraint functions can 
be written in the forms 

*l> » max max |h I, r - 1,2, ...R 
r t * 

L U (2) 

c k - c k - v k " 1 » 2 » •••* 

where h r 's are the R response functions, and C^, C^' are the prescribed 
upper and lower bounds of the kth constraint. We will restrict our con- 
sideration to cases wherein h f and are linear functions of the response 
(oi state) variables. 

The first step in using the indirect synthesis method is to calculate 
the limiting performance. Those portions of the syste-’ to be designed 
are replaced by generic (or control) forces {u ( t ) } , so that the equations 
of motion now become 

[M] {*} + [C) (x) + [K] {x} + IV] (u) - If] ( f ) (3) 

where [ K ] is the new stiffness matrix, which omits any stiffness con- 
tributions from the design elements, and [V] is the N x J coefficient 

matrix associated with the generic forces lu), J being the number of 
controllers. Such a general dynamical system is shown schematically 
in Fig. 1. Two types of elements are considered; structural elements and 
isolator elements. There may be any number of each; in particular M 
structural elements and J isolator elements could be Interconnected in 
an arbitrary fashion. A structural clement may represent a discrete 
mass point, a rigid body of distributed mass, or a flexible structure 
such as a framework or a shell. The isolator elements, similarly, can 
represent e'ther simple mechanisms without mass or models of more complicated 



devices. In general, the structural elementB constitute the prescribed 
portions of the system (l.e., the base structure and the elements to be 
isolated) and the Isolator elements are to be chosen in accordance 
with the design objectives. Note also that the control forces (u(t)}, 
being considered as explicit functions of time, do not contribute to the 
nonlinearity of the equations of motion even though they may have replaced 
nonlinear portions of the system, e.g. a nonlinear spring connection. 

Once Eqs. (3) are established, and the performance index v and the 
constraints are placed in the form of F.qs. (2), then existing limiting 
performance capabilities can be used to obtain the limiting performance 
characteristics which are subsequently employed in selecting the design 
parameters. 

As mentioned in Keference 1, the problem of establishing the 
coefficient matrices ( M ) , [ C J , [K] , (\), and ( F] in Eq. (3) 
for a particular dynamic system can be avoided by employing a general pi r- 
pose dynamic program to generate impulse responses at Isolator attach- 
ment points and then constructing h r , with Duhamel (convolution) in- 
tegrals. That is, since the system under consideration is linear, the 


response of the system excited by arbitrary forces u 


J 


can be obtained 


by superimposing responses to unit impulses placed sequentially at each 


of the isolator attachment points, i.e. 


J t 

h (t) - h (t) + l / g .( t - t)u (r)di 
r ru j-1 0 rj J 

(A) 

C k <t) - C k0 (O + - T) Uj (,)dT 


where g and g kJ 


are the appropriate system responses to a unit impulse 
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at the* <it tachment point ot tlu- Jih isolator, and h^ft), f^ft) nre 
the responses to lie I. input* f (t). The g and R. , must he generated 
by the Blructur.il dynamics program. Tlu* advantage of this approach lies 
in the fai t that .in existing general-purpose structural dynamic* .-ode 
will he used to perform ail of the system dynamics required for 
the limiting performance problem. However, in this work we are con- 
cerned with the situation whereby the system dynamics are to he solved by 
a limiting performance capability and rot by the general purpose progifiir. 
The general purpose program will be used to generate the coefficient 
mat r ices of Kqs . ( 1) . 


To illustrate the problems involved in preparing the equations for 
a limiting performance solution, consider the simple example of a three 
degrees of freedom spring-mass system as shown in Fig. 2a. Let it be 
desired to select the spring constant k 0 such that the maximum accelera- 
tion transmitted to am mass is minimzed while the three rattlespaces 
satisfy prescribed constraints. The k^, k, spring rates remain at 
their prescribed values. For ‘his example, the response functions 
that make up the performance Index are 
h ) - z v h 2 - z 2 . h 3 - z 3 

and the constraint filiations arc 


(: 1 " Z 1 - V ( 2 ‘ z 2 - *y ( 3 ‘ Z 3 ' f 
Therefore, Lqs. (3) become 


or 


tl' » max maxlz I, r ■ 1,2,3 
r t r 

ij/ ■ max (max | z. | , max|z 0 | , max|z-|] 
t 1 t i t 3 


_L . , ,.U 

C 1 - l ? 'l " z 2> - C 1 
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C 2 - '*2 *3 ^ - C 2 


c i : 1*3 - f l : c 3 

where C*‘, C*j , ...C^ .ire prescribed constants, i.e. prescribed constraint bounds. 
The equations of motion for tilts system nre easily found to be 


m^ 0 0 


m, 0 
0 m , 



0 


■kj k -k 2 

0 -k^ k 0 +k^ 



/■ > 


— 


X 1 


m l 

< 

X 2 


m 2 


X 3 

L * J 


1 

9 


(5) 


where x. 


- f, 1 ■ 1,2,3. Now replace the middle spring by a generic 


force u ? (Fig. 2b). This results in a new set of equations 


m l 

0 

0 


nip 0 



-k. 


-k. 


0 

0 

k„ 





— 


— — 


x i 


0 


“l 

< 

X 2 


1 

u 2 - - 

m 2 


X 3 


-1 


U> 

\ 


( 6 ) 


Equations ( 6 ) can be obtained either by applying the conditions of equilibrium 
to a conf igu/ation taken from Fig. 2b or by directly making the substi- 
tution U 2 ■ k^ (*2 ~ * 3 ) “ k^fx,, " * 3 ) ln eac *' r,,w k q • (5), which amounts 
to setting k,, * 0 and taking the coefficients of k., in row 2 of the [ K ] 
matrix as the column of the (V] matrix. Obviously, the first approach 
becomes impractical when the system is complicated or when the system is 
being analyzed by FE Codes since it contains unknown elements. In the 
following section, the second approach will be generalized and discussed 
further. 


DEVELOPMENT OF THE EQUATIONS OF MOTION 

As mentioned, in order to use most existing limiting performance 
compute. - capabilities, we must express the FE system equations of motion 
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In the form of Eq. (3) rather than In the usual FE form of 
F.q. (I). Standard FK codes will always generate the mass 


matrix ( M ) , the damping matrix |C) , and the coefficient matrix (F) which 
can be fed into the limiting performance codes with little or no modification. 
Thus it remains to create the new stiffness matrix ( K] and the cont r rller 
matrix [V) . Here we describe two technic s. 

Technique 1 

The replacement of an isolator to element k^ by a control force 
amounts to substituting in the equations of motion the relations 

Uj ■ k^ZajXj (7) 

where the summation is over the degrees of freedom that are connected 
with isolator and the a^'s are the klnematlcal factors that are 
associated with each degree of freedom. This lmrlies that Isolator k^ 
no longer contributes to the assembly stiffness matrix and Instead is re- 
placed by an additional matrix (Vj, whose individual columns consist of the 
coefficients of the term k^'a^x in the equations of motion. With the 
aid of the element stiffness ratrix v/ritten in assembly (global) 
coordinate , we can easily obtain the new stiffness matrix [K] and the 
controller matrix [V] from [K] . From the node numbering scheme of a 
FE code, we know which elements of [K] are contributed by Isolator k^ 
so that the new stiffness matrix [K] is obtained by simply removing 
those contributions. In the FE program, set the spring elements to 
be optimized temporarily equal to zero and the desired [K] is produced 


by printing the stiffness matrix. The controller matrix [V] is formed 


by taking as its column the coefficients of k^ in a row (or column) of 
[K] , the row (or column) number being that of a degree of freedom 
which is enacted by the isolator. 
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technique 


Rewrite Kq. (1) as 

iMHxt ♦ |i:]U» ♦ (K) ( x) - ( 1) < f ) - (V) { u > 


V V V 

11 12 * * 1J 

(F){f) - V 21 V 22 ... V 2J 


V V 

'..1 ••• >1 


N 1 N2 



( F J { f } -/ V 21 



t 

CM 

V 


r > 

v u 

c 

1 

1 

• 

> u 2 - 

V 2J 

• 


• 


• 

V- 


The [Kj is found as In technique 1. The columns of |VJ are obtained by 
noting that the elements of |V| are simply tin* Influence coelf Iclents 
of t ho "loads" (ui. Thus to obtain the influence coefficients of a 
control force u ( (i-e. to obtain the 1th column of |VJ), set all u's 
equal to zero. save u^ which Is set equal to 1 (fig. 3). The resulting 
load vector generator! bv FK program is the desired eclumn. In practice, 
one must he i ireful to properly account for the constraints when gen- 
erating the |K ! . jV| matrices t rom a FK code, l.e. It Is necessary to 
ascertain whether the matrices are printed before or after the dis- 
placement constraints have been applied to the problem. 


Three Degree of Freedom System Kxample 

Consider the problem of Fig. 2 which was studied above. To 
treat a particular case, let nij ■ 0.5, m 9 • 1. , tn^ * 2.5. ■ 10, k 2 ■ 20, 

k 3 - 30. Substitution of these values in Eq. (5) gives 
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In terms of global coordinates the •tlfiness |K ,] of the middle spring 
can he written as 



( 10 ) 


It follows from Eq. (10) that the stiffness of the middle spring contributes 
to K^f K , ^ , K^ ;t and K^. Removing these contributions from { K ) of Eq. (9) 
lesults In 



(ID 


and from the coefficients of elements in row 2 of (K., ] we find 



which Implies that u., ■ - z.^) » k n (x., - x ( ). 

To generate the desired matrices [M], (i<) , and [V] with a general 
purpose progtam using the second technique, first compute [M] as though 
a dynamics problem were to be run. Then run a static case to compute [Kj 
and [V] applying unit loads as shown in Fig. 4. 

PERFORMANCE INDEX AND CONSTRAINTS 

Now that the equations of motion are in the proper form for a limit- 
ing performance study, we proceed to establish the performance index and 
constraints . 
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Any linear combination ol ui > flfr.it li>m< , velocities, displacements, 
or control forces can be used a*. .1 performance index or constraint 
function, l.e. 

(h) or (C ) - IT] <x) 4 |Q) fx) + | H ) { x ) 4 ISHu) 4 IT) { f ) (12) 

where (P) , (Q) , (Kj , IS) and (" ) are prescribed coefficient matrices. 

Hie coefficient matrices In Eq. (12) t .in be easily formed when 
using a FE code. For example, consider the discretized system shown 
In Fig. 5. At each node, six degrees of freedom (DOF)f'i translational and 
3 rotational) are specified. Suppose in Fig. 5 that the accelerations 
along the y coordinate at node i and node j are to be constrained 
such that 

|ay t 4 by | < d (13) 

where a, b, d are prescribed constants. To place the constraint function 
a'y j 4 by j in • f orm of Eq. (12), we need only input the node numbers, 
the nodal i«)F, and the coefficients a, b, e.g. (1,2, a) and (J,2,b) where the 
2 stands for the y degrees of freedom. In most general purpose FE 
codes, there is a printable connect ivi tv array that identifies the nodal 
DOF to be constrained in terms ot the independent assembly DOF. Suppose 
nst ti,2) and (1,2) correspond to 1th and Ith assembly DOF. 

Then an appropriate kth row ol the [P| matrix in t.q. (12) will consist 
of the elements 


• a, j ■ b, k being the row number 

(14) 


P. - 0 for all f 4 I ,J. 
kt 

For this example the elements of the kth rows of the matrices [Q ) , (R), 
(S), |T] are all zero since the constraint function of Eq. (13) contains 
only accelerations. Similar reasoning permits Eq. (12) to be formed for 
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any acceptable type of performance Index and constraint function. 

Although we have chosen to describe a dynamic system by second 
order differential equations, Kqs. (3), some Uml.lng performance 
programs work with first order differential equations. The second 
order equations are convert ed Into a set of 2N first order linear 
differential equations 

{■} - (A] ( s ) + ( B ) { u } + [Dj { f } ( 15 ) 



is a state vector with 2N components, and the coefficient matrices 
(A), |B), [DJ are determined as functions of the coefficient matrices 
of Eq. (3). Also, Instead of Eq. (12), the performance Index and 
constraints are input to the limiting performance programs in the 
form 


(h) or (C) - [E] (s } + [Cl {u > + [H] {f } (16) 

where [El, [G] , and [Hj are prescribed coefficient matrices. Although 
the formation of Eq. (15) using Eq. (3) is a familiar manipulation, it 
is useful to cone ul nr in detail how Eq. (16) is obtained from Eq. (12). 
This is accomplished by solving Eq. (3) for (x) and placing the 
result in Eq. (12) to give 


{h} or {C} - ( [Q] - [P] [«-]){x) + ( [R] - [P] 

+ ([S] - [P] [~l~] ) (u } + ([T] + [P] [-£-])<f) (17) 
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A comparison of Fq. (IS) villi l.q. '!/) shows that 


in * fiQi H‘K-^-3 1 M - H-jHj 

|C] - [S] - ir][-~] 

and 

III] “ IT] * U'Jl "] 


( 18 ) 
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I'SF. OF SAP IV 


To implement the techniques described here the SAP IV computer pro- 
gram was selected. SAP IV (Reference ( ) is a popular general puroose program 
that can be altered locally with relative ease. The nropran required only 

slight modification in order to print out the desired matrices for Eqs. (3). 

In the following examples, all springs were input to SAP IV as elastic 
beam elements with appropriately chosen moduli and dimensions to secure 
the correct spring constants. The program allows for slave degrees of 
freedom to treat the rigid beams. It, however, ignores any loads 
applied at the slave nodes and thus these loads must be replaced by 
their static equivalents at the master node. 

Three Degree of Freedom System Example 

SAP IV was used to form the [M] , C K] , and (V] matrices for the 
three degree of freedom example of Fig. 2, with mass and spring values 
of the previous example. The results are given in Fig. 6. In its 
dynamic mode of opeution, SAP IV computed (M] . Then, in a static 
mode (K] and [V] were found by applying the unit loads of Fig. 4. * 

Five Degree of Freedom System Example 

As a more complicated example consider the automobile model of 
Fig. 7, which is considered in Reference 5. T he equations of rot ion for 
the configuration shown are 

Vl + k l (z l“ f ) “ k 2* z 2 ~ z l ■ J m 0 


" 2*2 + k 2 (z 2 " 2 1 " f e) - k 3 <*3 

19 ' k 2 (z 2 " Z 1 " I fl) f + V z 2 


z 2 ) + k 4 (s 


*4 + f 0) 


2 

21 

3 



- 0 


0 


4 




5 



If k., is replaced by u., so tli.it 


"2 * k 2 (l 2 - *1 ‘ 3 "> 


Chen we liave 


(KJ - 



0 

0 

0 

■ 

0 


-1 

0 

k 3‘ fk 4 

*L k 

3 4 

“ k 3 

” k 4 


+1 

0 

— k 

3 *4 

41 2 
-ii— u 

9 4 

0 

- — k 
3 *4 

, [v] - 

-£ 

3 

0 

“ k 3 

0 

k 3 

0 


0 

0 

" k <i 

3 *4 

0 

If +lf 

4 5 


0 


, {u} • u, (22) 


If In addition, k ^ is replaced by a control force so that 



As a numerical example take ■ tn^ - 2, m^ ■ 8, I ■ 0.5, ■ 1, k^ ■ k^ ■ 1, 


k2 ■ k^ ■ 2, k^ ■ 4, f. - 12, so that 
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In using SAP IV to generate these matrices, we first ran it in a 
dynamic mode to obtain [M] and [K] (displayed in Fig. 8) and then ran 
it as a static case by applying unit loads as shown in Fig. 9. This 
latter run yields the desired [K] and [V] matrices (Fig. 10). 

SUMMARY 

The techniques have been established whereby a general purpose 
structural analysis program can be used to form the equations required 
for available limiting performance capabilities (references 3,4,5). These 
techniques would permit a limiting performance study to be conducted on an 
arbitrary structural or mechanical system, which normally would involve 
finite element modeling for an analysis. The SAP IV general purpose 
computer program was employed to carry out a number of numerical 
examples to verify the formulation presented here. 
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Fig. 3 - Replacement of Generic force 1 > v I'nit Loads 
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Fig. 4 - Replacement of Generic Force by Unit Loads 
for System of Fig. 2 
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Fig. 5- Configuration for Forming Conatralnta and 
Performance Index 


23 



MASS MATRIX (DIAGONAL ELEMENTS ONLY) 
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Fig. 6 - SAP IV Generated Matrices for System of Fig 
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Fig. 8 - Matrices for Model of Fig. 



Fig. 9 Unit Force* for Model of Fig. 7 
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MASS MATRIX' DIAGONAL ELEMENTS ONLY) 
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